Cosmic Ray Positrons from Cosmic Strings 
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We study the spectrum of cosmic ray positrons produced by a scaling distribution of non- 
superconducting cosmic strings. In this scenario, the positrons are produced from the jets which 
form from the cosmic string cusp annihilation process. The spectral shape is a robust feature of 
our scenario, and is in good agreement with the results from the recent PAMELA and ATIC exper- 
iments. In particular, the model predicts a sharp upper cutoff in the spectrum, and a flux which 
rises as the upper cutoff is approached. The energy at which the flux peaks is determined by the 
initial jet energy. The amplitude of the flux can be adjusted by changing the cosmic string tension 
and also depends on the cusp annihilation efficiency. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Recent results from the PAMELA JL,] and ATIC 01 ex- 
periments have indicated an excess power of the cosmic 
ray positron flux compared to what is predicted from as- 
trophysical backgrounds alone. The power of the flux is 
observed to rise towards an upper cutoff which is in the 
range of 600GeV and falls off quite sharply above this cut- 
off. The data is in conflict with what is expected from 
astrophysical backgrounds. A peak in the spectrum of 
cosmic ray positrons in the energy range of 80 — SOOGcV 
had been suggested as a signature of dark matter anni- 
hilation in the galaxy a long time ago (the positrons are 
either produced from jets formed by the non-leptonic de- 
cay products J3j0, Q or by direct decay into electrons and 
positrons - for reviews see e.g. However, 

the specific form of the spectrum obtained from observa- 
tions is hard to reconcile with the predictions from simple 
dark matter annihilation models [11], although modified 
models have been proposed which are in better agreement 
with the data [T^. It is also possible that the positrons 
are due to nearby pulsars However, in the case of 

pulsars the decline of the flux at energies larger than the 
peak energy is not abrupt (see, however, [TJ] for a differ- 
ent point of view on this issue). 

In this Letter we would like to propose an alternative 
explanation for the positron excess which has nothing 
to do with dark matter. We investigate the possibility 
that the observed positron flux is due to jets from cusp 
annihilation of cosmic string loops. We are considering 
non-superconducting cosmic strings. It has been known 
for a long time that particle emission [T^ from cusps 
of cosmic strings leads to a spectrum which rises as a 
function of energy up to a cutoff set by the parameters in 
the cosmic string model and falls off quite sharply above 
this cutoff [16]. 

Since cosmic strings arise in a large class of particle 
physics beyond the Standard Model, our mechanism pro- 
vides a way to test physics beyond the Standard Model 



independent of the existence of low-energy supersymme- 
try. 

In the following section we review the basics of cosmic 
string dynamics and particle emission from cosmic strings 
which are important to understand our scenario. Section 
3 contains the computation of the cosmic ray positron 
flux. In Section 4 we give a brief discussion of the change 
in the positron spectrum during propagation through the 
galaxy. We conclude with a discussion of our results. 



II. PARTICLE EMISSION FROM A COSMIC 
STRING LOOP 

In this section we review the basics of cosmic string dy- 
namics and discuss the cusp annihilation mechanism by 
which cosmic string loops can emit high energy particles. 

Cosmic strings are one-dimensional topological defects 
which form during symmetry breaking phase transitions 
in a wide class of gauge theory models (see [13, [H, [l^ 
for reviews). If the gauge symmetry group at high tem- 
peratures is G, and the unbroken subgroup below the 
transition temperature is H, then the criterion for the 
existence of cosmic strings in the theory is 



Hi(A^) ^ X, 



(1) 



where Ad is the vacuum manifold of the theory below 
the transition temperature and X is the trivial group. If 
the group G is simply connected, then M. = G/H. This 
criterion is satisfied in a large class of particle physics 
theories beyond the Standard Model. 

In particle physics theories admitting the existence of 
cosmic strings, such strings inevitably arise during the 
symmetry breaking phase transition ^20]. By causality, 
the point in M. which the order parameter describing the 
phase transition takes at temperatures lower than the 
transition temperature must be uncorrelated on length 
scales larger than the Hubble radius H^^ ^ where H is 
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the cosmological expansion rate. Hence, there is a prob- 
abihty of order 1 that one string wiU traverse any partic- 
ular Hubble volume after the phase transition. Typically 
(in particular if matter above the transition temperature 
is in thermal equilibrium), the initial separation of the 
cosmic strings will be microscopic. The above causality 
argument applies at all times subsequent to the transi- 
tion time. Thus, in a theory which admits cosmic strings, 
then at all times after the phase transition (in particular 
at recent cosmological times) a network of cosmic strings 
with separation no greater than the Hubble radius will 
be present. 

Cosmic strings arising in gauge field theories must be 
closed, i.e. either string loops of "infinite" strings (de- 
fined as truly infinite strings or string loops with cur- 
vature radius larger than the Hubble radius. Thus, the 
system of cosmic strings at any time after the phase tran- 
sition will consist of a network of "infinite" strings and 
a distribution of string loops. Both analytical arguments 
detailed i n [TtI fisl . and detailed numerical simula- 
tions [U I22I |23|| have shown that the network of infinite 
strings approaches a "scaling" solution characterized by 
a string correlation length which is a fixed fraction of the 
Hubble radius at all late times. Roughly speaking, we 
can view the long string network as a random walk with 
step length of the order of the Hubble radius. 

The scaling solution implies that the total length in 
the long strings is decreasing. This decrease is realized 
by the inter-commutation of long strings. Such inter- 
commutations produce string loops. Thus, at any time 
t, there will be distribution of loops. Neither numerical 
simulations nor analytical studies at this point agree on 
the exact nature of the loop distribution (see e.g. Q 
for recent progress). We will be using a simple one-scale 
model for the distribution of loops which is based on the 
assumption that loops at time t form at a fixed fraction 
of the Hubble radius. 

Once formed, string loops decay predominantly by 
gravitational radiation. The rate of gravitational radi- 
ation is governed by the string tension /i, namely 

7?-7Gm, (2) 

where G is Newton's gravitational constant and 7 is a 
numerical constant whose value is of the order 10^ psj . 

Taking into account the redshift of the number density 
of string as well as the formation scenario and decay rate 
of cosmic strings discussed above, we obtain the following 
distribution of string loops at time t (see e.g. [2^]): 

n{R,t) = KR-^r^ , R > teq (3) 
n(i?,t) = KR-'^'/^rhl/q^ , -1G^it < R < teq 

n{R,t) = K{-fGn)-^/^t-hll^ , R < -iG^it, 

where k is another numerical constant which depends on 
the details of the cosmic string scaling solution, and t^q is 
the time of equal matter and radiation. In the above, we 
are assuming 7G/i < teg- The first line in represents 



loops which were produced after the time of equal matter 
and radiation, the second line loops generated before t^q 
which will survive gravitational radiation for more than 
a Hubble expansion time. The last line represents loops 
which are in the final stages of decay by gravitational 
radiation. 

From the distribution ([3]) it can be seen that the energy 
density in string loops is dominated by loops of radius 
about R ^ '^Gfit. These loops will, as we show below, 
also dominate the positron emission from strings. 

Cosmic strings as two-dimensional world sheets 
x'^ (cr, r) (where r is a world sheet time coordinate and 
a labels the spatial world sheet coordinate) are solutions 
of the Nambu-Goto equations, the same equations which 
describe fundamental strings. Since cosmic strings have 
relativistic tension, they will oscillate. It can be shown 
^2?| that "cusps" generically occur on string loops (at 
least once per oscillation time) . A cusp is a point on the 
string where 

x^ = 0, (4) 

where a prime indicates the derivative with respect to cr. 

Geometrically, a cusp corresponds to a "spike" on the 
string (see Figure 1). Since a cosmic string has a finite 
width w whose magnitude is of the order of at a 

cusp the two segments of the string at either side of the 
cusp overlap. By expanding the solutions of the Nambu- 
Goto equations about a cusp it can be shown [28i] that 
the length of the overlap region is 

I, ^ w^/^R^/' , (5) 

where R is the radius of the string loop. 




W 



radiation jets 

FIG. 1: Sketch of a cusp on a cosmic string loop. The mean 
curvature radius of the string loop is R, the string width is 
w, and the length of the overlap region is Ic- 

There is nothing that prevents this overlap region from 
annihilating into particle excitations of the Higgs and 
gauge fields which the string is made up of. Assuming 
that the entire overlap region decays, one gets a power of 
particle radiation given by (see also [2§|) 

Pc - filcR'^ /xwi/3i?-i/3 . (6) 

This particle decay rate is - for string loops of macro- 
physical radius R - very small compared to the power 
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radiated into gravitational radiation which is 



(7) 



(see <^). Nevertheless, as was discussed in [1^ (see also 
[3011 )■ cusp annihilation can contribute a significant frac- 
tion to the cosmic ray flux, assuming that a substantial 
fraction of the string overlap region decays [i^ 

The expected ultra-high energy neutrino flux from cos- 
mic strings with a tension given by the scale of Grand 
Unification was studied in [igI . [soj , and the correspond- 
ing 7-ray signatures were analyzed in [s^ . These works 
were performed under the assumption that the distribu- 
tion of cosmic string loops scales as described by 
The analysis was extended in [3^ to the case of a non- 
scaling loop distribution (numerical evidence for such a 
non-scaling loop distribution came from the analysis of 
[33|). However, most studies of cosmic string dynamics 
favor a scaling distribution of string loops (see e.g. [24|). 
and thus in this paper we will assume such a scaling dis- 
tribution. 



III. POSITRON FLUX FROM A SCALING 
NETWORK OF COSMIC STRINGS 

As discussed above, the power of energy loss from a 
cosmic string loop of radius R is 



Pr 



5/6^-1/3 



(8) 



where e is the efficiency factor discussed above and we 
have made use of the fact that the string width w is pro- 
portional to the inverse square root of the string tension 
A*- 

The primary decay products from cusp annihilation are 
quanta of the scalar and gauge fields which make up the 
cosmic strings. These quanta, in turn, will decay into rel- 
ativistically moving standard model particles which will 
form jets. Following the discussion in [l^, we take the 
primary energy of a single jet to be m/. In this case, the 
number of jets N formed per unit time is 



N 



en'^/^mJ^R-^/^ = PR 



-1/3 



(9) 



where the last step defines the quantity P. 

A single jet leads to the following spectrum of energies 
(numberper energy interval) of neutrinos resulting from 
the jet |35| : 



dN 15 _wll ^ , _i/2 2 _ 

dE 16 f ^ 3 3 



3/2\ 



where 



E 
ruf 



< 1. 



, (10) 



(11) 



We will take the same formula to give the energy spec- 
trum of all stable leptons resulting from the decay, in 
particular the spectrum of positrons. 



Particle physics models admitting non- 
superconducting cosmic string solutions will have 
string cusps decaying at all times. Neutrinos produced 
from cusp decay will travel cosmological distances (see 
e.g. [lB|)j whereas positrons are absorbed and lose their 
energy on super-galactic scales. To obtain the neutrino 
flux at energy E, we have to integrate over all times 
t the flux of particles emitted at time t with energy 
E{z{t) + 1). For electrons and positrons, we only have to 
integrate over times which are smaller than the current 
time to by less than the "containment time" of electrons 
and positrons, respectively, in the galaxy. According to 
0, the containment time of positrons is of the order 10^ 
yrs,, i.e. much longer than the time it would take light 
to travel through the galaxy. The containment time 
corresponds to a redshift of Zc ~ 10^"^. 

When computing the expected positron flux from cos- 
mic string cusp annihilations we must therefore impose 
several cutoffs. First of all, only positrons emitted at 
redshifts smaller than Zc will contribute. Secondly, only 
string loops located inside the galaxy may be considered. 

The general expression for the differential energy flux 
F{E) of cosmic ray positrons from cosmic string cusp 
annihilations is 



F{E) 



dt{z{t) + l)-^f{iz{t) + l)E,t), (12) 



where / is the differential flux per unit time of positrons 
emitted at time i, the "injection spectrum" . The redshift 
enters in two places. Firstly, the injection number density 
is redshifted, and, secondly, the energy of a given positron 
redshifts. 

The injection spectrum is obtained by integrating over 
all cosmic string loops present: 



dN . 



f{{z{t) + l)E,t) = (z(t) + l) — U^(,+i)B/,„^ 
dRn{R,t)PR-^^^ , 



(13) 



where the string number density n{R, t) and the constant 
P have been defined previously. The first factor on the 
right-hand side of this equation is the Jacobean factor 
obtained by transforming between final energy E and 
injection energy. 

As follows from recalling the loop distribution ([3]), 
the integral over R is dominated by loops of radius 
R ~ 'yG^t. If we consider the string scale to be at least 
a couple of orders of magnitude smaller than the scale of 
Grand Uniflcation, string loops of radius R ^ jGfito will 
still be present today, and their number density will be 
such that many string loops of such radius will be located 
within our galaxy. 

To estimate the amplitude and shape of the spectrum, 
we flrst insert the injection flux (|13p into ll2p and perform 
the integral over loop radii. The integral is dominated 
by the value R = jG^t. A good estimate of the result is 
obtained by integrating over the loops with radii in the 
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range ^G^t < R < t^g. The result is 
F{E) ^ P,.{^Gf,y/'tl/' 

2,-2 ,-11/6 



(14) 



dN. 



\x={z+l)E/mf 



Next is the integral over time which can be simplified by 
using the integration variable 



z = z{t) + 1 , 



(15) 



To obtain an estimate of the flux, we use for dN/dE 
the final term on the right-hand side of (|10p . After a 
couple of lines of algebra (and in particular plugging in 
the definition of P from®) we obtain 



E^F{E) 



11/6^-3/4 
eg 



(16) 



where the final factor comes from the range of integration 
over t. 

Inserting numbers into (I16|) and expressing the result 
in terms of the units which experimentalists use we get 



E^F{E) ~ ev-f-^/^{-fGfi)~^ 

™/ Gcv( 



(17) 



z,10""m~2gg(,-iQgv2 



The specific signature of our predicted cosmic ray 
positron flux is the power law increase oiE^F{E) ~ E^l^ 
and the sharp cutoff at an energy scale set by the ini- 
tial jet energy to/. In contrast, the background flux of 
positrons (multiplied by E^) from astrophysical sources 
is predicted to be slightly decreasing in the energy range 
between lOGeV and lOOOGeV. 

Let us first give a rough analytical treatment of the 
predicted positron to electron flux ratio. Both fluxes are 
a superposition of background and cosmic string-induced 
fluxes, and we will use the subscripts hg and cs, respec- 
tively, to denote these two contributions. The flux ratio 
$ is 



E^F{E)- 



E^F{E)- +E^F{E)- 



(18) 



Assuming that the electron flux is dominated by the 
background, we obtain 



^ = R, 



bg 



E^F{E); 
E^F{E) 



(19) 



bg 



where Rbg is the background flux ratio. Since E^F{E)'j^g 
is roughly constant, we see that the flux ratio in the cos- 
mic string model is predicted to be equal to the back- 
ground value at low energies and gradually shift to scal- 
ing as 



(20) 



at higher energies. At energies close to the cutoff value 
TO/, the spectrum again flattens out because terms in (jlOp 
scaling differently than the x~^/^ term which we focused 
on will become important. 

From the PAMELA data [H we can read off a slope 
which is rising to about 0.5 at energies between 50 and 
lOOGeV. From the ATIC data, a slope of close to 1 is 
inferred at energy scales between 300 and GOOGeV. 

To obtain a better idea of the fit of our model, we have 
evaluated the predicted positron flux numerically, keep- 
ing all of the terms in (fTUl) . Our results are plotted in 
Figures 2-4. The numerical results also include the pro- 
cessing of the spectrum during propagation as discussed 
in the following section. 



IV. PROPAGATION OF THE POSITRONS 

Positrons will lose energy not only because of red- 
shifting, but also because of interactions during their 
propagation from the source to us through the galaxy. 
We consider a standard diffusion model for the propaga- 
tion of positrons in the galaxy. 

To begin with, let us recall the physical processes 
which affect the propagation of charged particles in the 
galaxy. Firstly, when a charged particle travels through 
the galaxy, its movement can be affected by the galactic 
magnetic field. Although the magnetic gyro-radius of a 
particle is usually very small, this particle can still pos- 
sible to jump to near-bye field lines due to the tangled 
magnetic field and so could change its orbit. We usually 
model this process with a diffusion equation. Secondly, 
during the propagation of a positron, the particle loses 
energy because of inverse Compton and synchrotron pro- 
cesses. These two factors are the most imp ortant ones. 
For a detailed study, we refer to Ref. [30|. In the fol- 
lowing we focus on the above-mentioned two processes, 
especially the energy loss 

Neglected other effects which are present in addition 
to the two mentioned in the previous paragraph, and 
assuming a spherically symmetric diffusion process, we 
obtain the following propagation equation for the flux F 
of charged cosmic ray particles: 



d_ 



F ^ D{e)\/^F+^jL{e)F] +Q{e,x) , (21) 



where e is defined as a dimensionless energy variable e = 
with a the scale factor of the universe, D is the 



a GeV 

diffusion coefficient, L is the energy loss rate and Q is 
the source term. 

In models in which dark matter annihilation is the 
source of the positron excess, the production of positrons 
is dominated by the annihilation of dark matter particles 
today. Therefore, usually only steady state solutions of 
Eq. (|2ip (in which he left hand side of the equation van- 
ishes) are considered, as analyzed for example in Ref. [13] ■ 

However, in our model the source of positrons does not 
scale in time as the background density, and hence the 
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resulting flux will not be steady state-like. As discussed 
in the previous section, we need to integrate the flux 
equation from the earliest moment from which positrons 
will still reach us today. This time is the containment 
time of positrons in the galaxy which corresponds to a 
redshift ~ 10~^. 

In the following, we will neglect the diffusion term in 
(PTjl . In models with a string tension significantly smaller 
than that given by the scale of Grant Unification, the 
separation of strings is much smaller than the radius of 
the galaxy. Many of these strings are a distance away 
from us which is smaller than the diffusion radius. Hence, 
we argue that we can treat the flux as quasi-homogeneous 
and hence neglect the diffusion term. 

Thus, we consider the simplified propagation equation 



dt 



+ + !)-'/((«(!) + l)E,f) (22) 



To solve this equation, we apply a perturbative approach 
and separate the flux F into inflnitely many components. 



F = Fo + Fi 



(23) 



Each component satisfies its own propagation equation 
as follows, 

= 



and, more generally. 



(24) 



After solving these equations one by one, we then sum 
up all the component to obtain the result 

F(e) = f^dta^f+ r dt^L{€) f dt'a^f+ ....(25) 

J to J to Jig 

By changing the time integrals to all the entire integral 
range, we can obtain the factor 1/n! for the n-th compo- 
nent. Eventually, we obtain the following formal solution 
for the flux: 



Fie) = Fo(e,to) xexp ^^^^ 



/;;dtfL(e)Fo(6,t) 



, (26) 



where _Fo(ej^o) is exactly what we have studied in Sec. 
Ill without considering the energy losing effect. 

In a realistic model, we consider energy loss through 
synchrotron emission and inverse Compton scattering. 
As introduced in Rcf. [38], the process can be param- 
eterized as 



Lie) = - 

TE 



(27) 



with the energy-loss time te — 10 s. After combining 
Eq. (16) and the energy- loss parametrization, we can 
derive the factor in the exponential term of Eq. (pS]) . It 
is 



Zo+Zc 



, 3to / _3 

az [z ^ — z 

4te 



l)E- 



ZrE- 



-3/2 



9^0 ZcE 
"32TB IGeV 



(28) 



in the low energy regime. Correspondingly, we obtain the 
following approximate form of the flux 



FiE) ~ Fo(£;)exp{- 



9to z^E 
32te IGeV 



} 



(29) 



The exponential factor describes the energy losing of the 
positrons when they are passing through the galaxy. 

From the above result, we learn that because of the 
smallness of the containment time of positrons in the 
galaxy, the energy loss due to interactions is insignificant 
for low energy positrons. For higher energy positrons 
the energy loss becomes more important. This leads to 
a slight smoothing out of the delta function-like upper 
cutoff in the predicted flux. In addition, depending on the 
parameters, the energy corresponding to the maximum 
of the flux may be smaller than mf. 

To obtain a better idea of the fit of our model, we have 
evaluated the predicted positron flux numerically, keep- 
ing all of the terms in (10) and taking into account of 
the energy loss effects on the propagation of positrons 
discussed in this section. Our results are plotted in Fig- 
ures m [3] and 01 In the figures, we take three groups 
of parameters for the model as shown in the captions of 
these figures. We have chosen the parameters v — 13 
(determined by numerical simulations of cosmic string 
evolution [1^), the containment time Zc = lO^'^, and the 
energy-loss time te — lO^^s. 

Whereas the PAMELA experiment only shows a rise 
of the flux ratio as a function of energy, the ATIC exper- 
iment which probes the spectrum of positrons to higher 
energies shows a sharp upper cutoff at an energy of about 
600 GeV. Thus, matching the ATIC data leads us to pre- 
fer a higher value for the initial jet energy m/. 



V. CONCLUSIONS AND DISCUSSION 

In this Letter we have proposed a new explanation for 
the observed excess of positrons over electrons n the cos- 
mic ray flux at energies between 10 and 600GeV. In our 
model, the source of the cosmic ray positrons is radiation 
from cosmic string cusp emission. Our model does not 
assume that the particle physics model manifests super- 
symmetry with low-scale supersymmetry breaking. In 
contrast, it assumes the existence of linear topological 
defects. 

The spectral shape which we predict is insensitive to 
the details of the cusp annihilation process and is thus a 
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FIG. 2: The predicted positron flux from cosmic string cusp 
decay for various values of the string tension G/i and the initial 
jet mass of m/. The values of the other parameters were 
chosen Xohe v = 13 and efficiency factor e = 1. 
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FIG. 3: Comparison of the cosmic string model predictions for 
the flux ratio $ with the data from the PAMELA experiment, 
for the same choices of parameters as in the previous figure. 

robust prediction of our model. The position of the peak 
of the flux is determined by the initial jet energy m f . The 
amplitude of the spectrum, is not a robust prediction of 
our model. It depends sensitively on both the cosmic 
string tension and on the efficiency factor e of the cusp 
annihilation process. In our plots, we have fixed e = 1. 
From the analysis in Section III it follows immediately 
that the factor which determines the amplitude of the 
flux is e(7^^/^G'/i)~^. This is the factor which can be fixed 
from the recent positron flux observations, assuming that 
our mechanism is the source of the excess. 

From our analysis we can learn another important 
lesson: for flxed value of m/, a model with non- 
superconducting cosmic strings predicts a positron flux 
with a shape given by our analysis. Even if the observed 
flux is not due to strings, we get an upper bound on the 



FIG. 4: Predictions for the flux ratio $ at higher energies, for 
the same parameter values as in the previous two figures. 

quantity e(7^^/^G/i)~^. As our results show, for e ^ 1 
this is a bound which rules out many models with low 
energy scale strings. 

Let us add some more comments on the sensitivity of 
our predictions to the value of the efficiency factor e. If 
we were to use to value of Ic given by [3l[, which takes 
into account jjffects which were not included in the initial 
analysis of [15], and is 

I, ^ w^l^R^I^ , (30) 

then the predicted amplitude of the flux decreases by a 
factor of [wji^Y^^ ■ Moreover, back-reaction effects on 
cusp formation are still not included completely in [sij . 
and thus the actual amplitude may even be lower. The 
assumptions we make about the efficiency of cusp annihi- 
lation will change the value of G/i for which the amplitude 
of the positron flux agrees with observations. 

There are large classes of particle physics models be- 
yond the Standard Model which predict the existence 
of cosmic strings. Cosmic strings are also predicted in 
many inflationary universe models based on superstring 
theory [s^ (for reviews see e.g. [i^). A network of cos- 
mic strings will also remain in the string gas cosmology 
model i41i, '4^ (for a recent review see (43j). 

In order for our model to be consistent with the ab- 
sence of an excess in the cosmic ray anti-proton spectrum 
[i^ l we require the cosmic strings to decay predominantly 
leptonically. 

In light of the recent positron data, our work motivates 
a closer look at the mechanism of cusp annihilation. Any 
improvement in our understanding of this process would 
lead to a much improved predictive power of our analysis. 
Another issue which merits re-visiiting is the determina- 
tion of the initial jet mass rrij resulting from cusp anni- 
hilation. Another important issue is to determine which 
cosmic string models can lead to predominantly leptonic 
jets. On the experimental side, it is interesting to ex- 
plore ways to distinguish between the proposed scenarios 
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to explain the positron excess (see e.g. j45|). 
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